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Introduction 


The  question  was  raised  as  to  the  suitability  of  using  a  95% 
nitromethane  (NMe)  and  5%  diethylamine  (DEA)  liquid  explosive, 
currently  in  use  by  the  Royal  Armament  Research  And  Development 
Establishment  (RARDE),  as  a  quick  turn-around,  easy-to-load,  and  uniform 
explosive  loading  for  developmental  liners  under  the  Advanced  Chemical 
Energy  Warhead  Program.  There  is  precedent  for  using  substitute 
explosives  during  liner  materials  investigations  and  liner  design 
development,  primarily  to  ensure  consistent,  precise  loading  quality  and  as 
well  as  for  convenience  and  lower  cost.  The  use  of  alternate  explosives, 
however,  is  only  suitable  to  the  point  that  the  substitute  explosive  does  not 
significantly  alter  the  conditions  under  which  the  intended  liner  is  to  be 
employed,  or  adversely  affect  those  material  and  liner  parameters  under 
investigation. 

In  this  analysis,  the  NMe  liquid  explosive  will  be  substituted  into  a 
BRL  3.3  inch  42  degree  conical  copper  lined  charge.  The  original  charge  is 
Octol  loaded.  Since  the  NMe  is  less  energetic,  it  is  expected  to  give  a 
reduced  jet  tip  velocity,  less  jet  mass,  and  lower  overall  jet  kinetic  energy. 
This  would  naturally  follow  from  the  use  of  a  lower  energy  density 
explosive.  The  liquid  explosive  load  will  then  be  modified  by  a 
combination  of  liner  subcalibering  and  detonation  front  wave  shaping  to 
investigate  how  the  original  jet  characteristics  may  be  recovered. 
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Approach 


The  most  precise  approach  to  analyzing  the  performance  tradeoffs 
between  the  two  explosive  loads  is  to  use  a  hydrocode  simulation.  The 
hydrocode  allows  accurate  modeling  of  the  detonation  and  gas  pressures 
behind  the  liner  over  the  entire  collapse  duration.  The  effects  of  charge 
confinement  and  rarefaction  shock  waves  on  the  collapse  velocity  are  also 
readily  apparent  in  the  simulation.  Placing  tracer  elements  along  the  liner 
contour  allows  the  collapse  velocity  to  be  monitored  and  the  collapse  angle 
to  be  calculated.  Once  the  elements  reach  the  axis,  application  of  the 
hydrodynamic  equations  for  jet  formation  readily  follows. 

Without  the  use  of  a  hydrocode  to  determine  the  liner  collapse 
velocity,  as  in  this  analysis,  an  approximation  is  required.  Two  possible 
formulations  present  themselves.  The  first  is  to  assume  a  pressure  history 
behind  the  liner,  beginning  with  a  peak  at  the  detonation  pressure,  and 
then  diminishing  linearly  or  quadratically  to  some  final  value,  based  on  the 
charge  confinement.  Not  knowing  with  confidence  how  confinement  and 
rarefaction  waves  affect  the  pressure  history,  a  second  approach  was 
formulated  based  on  the  Gurney  equations.  The  Gurney  equations  provide 
a  final  plate  velocity  which  accounts  for  both  charge  height  and  charge 
confinement.  The  first  problem  with  this  approach,  however,  is  that 
Gurney  velocities  are  based  on  experiments  with  parallel  plates  and  only 
provide  final  velocities  for  the  entire  plate,  not  discrete  sections. 

Application  of  this  to  axi-symmetric  conical  shaped  charge  liners  clearly 
violates  this  condition.  This  approach  also  does  not  provide  an  acceleration 
history  of  the  liner  element,  and  it  is  well  known  that  the  apex  of  the  liner 
generally  reaches  the  axis  prior  to  realizing  its  full  velocity  potential.  In 
addition,  the  angle  of  incidence  of  the  explosive  detonation  front  affects 
the  acceleration  time,  although  it  generally  does  not  affect  the  final 
velocity  of  the  plates. 

It  is  true  that  the  Gurney  velocities  are  experimentally  determined 
on  a  geometry  radically  different  from  a  shaped  charge.  However, 
intuitively,  the  shaped  charge  application  of  the  Gurney  equations  must  be 
a  more  general  case  of  the  actual  experiments.  In  other  words,  they  are 
related. 


The  Model 


The  acceleration  of  the  liner  elements  is  based  on  the  Gurney 
equations  for  an  asymmetric  explosive  sandwich.  The  liner,  explosive 
charge,  and  confinement  are  divided  into  100  discrete  segments,  which 
facilitates  the  calculation  of  charge  to  mass  and  charge  to  confinement 
ratios.  Working  in  terms  of  charge,  liner,  and  confinement  masses  neglects 
any  effects  of  the  cylindrical  geometry  in  the  shaped  charge.  Figure  1 
presents  the  applicable  equations  as  taken  from  Reference  1. 


ASYMMETRIC  CONFIGURATIONS 
Open-faced  Sandwich: 


Asymmetric  Sandwich: 


Symbols : 


M 

C 

v 


Gurney  energy,  E 


total  metal  mass  _  N  total  tamper  mass 

total  explosiv-  mass  '  C  =  total  explosive  mass 

metal  velocity 

kinetic  energy/unit  explosive  mass 


Figure  1  Gurney  Equations 
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Each  liner  element  is  assumed  to  achieve  85-90  %  of  its  full  Gurney 
velocity  within  the  time  it  takes  the  release  wave  to  arrive  radially  from 
the  charge  confinement.  The  release  wave  speed  is  assumed  to  be  the 
explosive  detonation  speed.  This  formulation  allows  for  the  condition 
where  liner  elements  reach  the  axis  prior  to  full  acceleration.  The 
acceleration  history  is  assumed  to  be  exponential,  beginning  with  zero 
velocity  when  the  shock  front  arrives  at  the  element  rising  to  85-90%  of 
maximum  velocity  by  the  time  the  release  wave  reaches  that  element  from 
the  charge  perimeter,  and  then  100%  Gurney  velocity  is  realized  at  an 
effectively  infinite  time.  The  85-90  %  Gurney  velocity  term  in  the 
formulation  was  chosen  because  it  fit  jet  tip  velocity  data  for  known 
shaped  charges.  For  the  octol  loaded  3.3  inch  charge  in  this  analysis,  using 
90%  gives  the  best  results. 

The  model  also  incorporates  plane  wave  detonation  fronts,  as  well  as 
point  and  ring  initiation  at  a  defined  charge  head  height  behind  the  apex  of 
the  liner.  When  a  liner  is  modeled  using  a  rounded  apex,  it  is  best  to  cut 
the  apex  off  and  ignore  it.  In  this  region  the  accuracy  of  charge  to  mass 
ratio  calculations  are  questionable,  the  liner  material  has  very  little  run 
distance  to  the  axis  so  it  usually  has  a  strong  reverse  gradient,  resulting 
only  in  a  massive  jet  tip.  This  model  will  not  accurately  handle  these 
rounded  apexes,  which  behave  more  truly  like  hemispherical  liners. 


Results 


Figure  2  presents  the  3.3  inch  BRL  standard  shaped  charge  as  taken 
from  reference  2,  which  was  the  baseline  charge  for  this  analysis.  Figure  3 
shows  the  computer  mesh  with  100  liner,  explosive,  and  confinement 
elements.  The  charge  is  Octol  loaded,  which  has  a  Gurney  velocity  of  2.83 
mm/usec,  detonation  velocity  of  8.5  mm/usec,  and  a  density  of  1.82  g/cc. 
Point  initiating  this  charge  at  a  head  height  of  67.6  mm  gives  the  following 
jet  parameters  according  to  the  model: 


Baseline  Octol  Charge  (BOC) 


Vtip  (element  13);  mass 
V  (element  100);  mass 
Jet  mass 
Jet  energy 


8.31  mm/usec  ;  .177  gr 
0.33  ;  2.77 
71.6  g 

.543  megajoules 
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Figure  3.  Computer  Mesh  of  Baseline  Octol  Charge 
(Only  one  half  of  axi-symmetr ic  shape  shown) 


The  jet  tip  has  a  slight  reverse  gradient  from  element  1  back  to 
element  13,  where  the  momentum  averages  out.  Actual  test  data  for  this 
charge  gives  a  tip  velocity  of  8.30  mm/usec. 

The  NMe  liquid  explosive  has  a  detonation  velocity  of  6.3  mm/usec, 
and  density  of  1.1  g/cc.  The  Gurney  velocity  for  this  explosive  is  unknown 
to  us,  so  it  was  estimated  based  on  a  linear  extrapolation  from  Octol 
through  TNT,  which  is  closer  in  properties  to  the  NMe  than  Octol  is.  The 
Gurney  velocity  is  assumed  to  be  a  function  of  the  pressures  behind  the 
liner  as  it  accelerates.  These  pressures  should  be  a  function  of  the 
explosive  density  times  the  detonation  velocity  squared.  Knowing  these 
parameters  for  Octol  and  TNT  (V  det  6.95,  density  1.65,  V  Gurney  2.37) 
yields  an  estimate  of  2.04  for  the  Gurney  velocity  for  the  NMe.  Running 
the  same  charge  with  NMe  yields  the  following  results: 

Baseline  Liquid  Charge  (BLC) 

Vtip  (element  5);  mass  5.67  mm/usec  ;  .128  g 

V  (element  100);  mass  0.19  ;  2.76 

Jet  mass  69.3  g 

Jet  energy  .206  megajoules 

There  has  been  a  significant  drop  in  jet  tip  velocity  and  jet  energy. 

About  10%  of  jet  mass  is  lost  in  the  forward  portion  of  the  jet,  but  recovers 
somewhat  in  the  rear  portion.  The  overall  jet  mass  is  nearly  the  same. 
However,  useful  jet  mass  (above  2  km/sec)  reduces  from  element  84  in  the 
baseline  charge  to  element  72  in  this  charge. 

Since  the  jet  tip  velocity  dropped  so  dramatically,  no  run  was  made 
using  a  plane  wave  detonation  front  to  simulate  more  explosive  head 
height.  Something  more  drastic  was  needed  to  bring  jet  tip  velocity  up. 
Going  to  ring  initiation  allows  the  detonation  front  to  sweep  by  the  liner 
faster,  artificially  increasing  the  detonation  velocity,  resulting  in  a  lower 
collapse  angle.  This  should  increase  tip  speed.  However,  it  wiil  be  at  a  cost 
of  lower  jet  tip  mass. 

Initiating  the  liquid  charge  at  a  head  height  of  20  mm  and  40  mm  off 
axis,  or  ring  initiation  at  the  charge  perimeter,  gives  the  following  results: 
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Ring  Initiated  Liquid  Charge  (RILC) 


Vtip  (element  1);  mass 
V  (element  100);  mass 
Jet  mass 
Jet  energy 


8.39  mni/usec  ;  .069  g 
0.19  ;  2.75 
66.7  g 

.216  megajoules 


Jet  tip  velocity  recovers  to  that  of  the  baseline  charge,  but  the  tip 
mass  is  about  60%  of  the  point  initiated  liquid  charge.  Jet  mass  recovers 
back  at  element  39,  but  then  begins  to  drop  off  by  about  10%  back  to 
element  73.  Total  jet  energy  is  still  less  than  half  of  the  baseline  charge. 

The  liquid  charge  was  then  subcalibered  as  shown  in  Figure  4.  The 
charge  diameter  is  now  140mm  rather  than  83  as  in  the  baseline  charge. 
Confinement  in  this  charge  is  5  mm  thick.  This  charge  was  ring  initiated  at 
a  charge  height  of  57  mm  and  a  radius  of  70  mm. 

Subcaliber  Ring  Initiated  Liquid  Charge  (SRILC) 

Vtip  (element  1);  mass  8.34  mm/usec  ;  .070  g 

V  (element  100);  mass  4.66  ;  .653 

Jet  mass  32.1  g 

Jet  energy  .616  megajoules 

Jet  tip  velocity  is  maintained  and  jet  energy  increases.  However,  jet 
mass  is  cut  in  half  and  the  tail  velocity  is  too  high  for  practical  application. 
This  charge  was  then  modified  by  tapering  the  confinement  down  at  the 
base,  in  hopes  of  bringing  the  jet  energy  and  tail  velocity  down.  Figure  5 
shows  this  geometry.  The  charge  diameter  at  the  apex  remains  140mm 
and  tapers  to  120  mm  at  the  liner  base.  The  same  initiation  is  used  as  last 
time. 


Subcaliber  Ring  Initiated  Tapered  Liquid  Charge  (SRITLC) 


Vtip  (element  1  );  mass 
V  (element  100);  mass 
Jet  mass. 

Jet  energy 


8.45  mm/usec  ;  .068  g 
3.21  ;  .911 
36.3  g 

.542  mcgajoules 


The  jet  parameters  are  moving  in  the  right  direction,  but  are  clearly 
not  going  to  match  the  octol  load.  It  appears  that  any  more  tapering  will 
only  reduce  the  jet  energy  below  the  baseline  charge,  despite  the  still 
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Figure  4.  Subjalibet  Liquid  Charge 
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Figure  5.  Sutma. i.  iber  Liquid  Charge 


elevated  tail  velocity.  Tip  velocity  is  good,  yet  total  jet  mass  is  still  half  of 
the  baseline  charge. 

To  see  how  more  tapering  affects  the  jet  parameters,  the  charge  was 
given  the  most  extreme  taper  possible  —  140  mm  at  apex  and  83  mm  at 
the  base,  as  shown  in  Figure  6.  Initiation  is  at  a  radius  of  70  mm  at  a 
head  height  of  65  mm.  The  results  are: 

Subcaliber  Ring  Initiated  Extremely  Tapered  Liquid  Charge  (SRIETLC) 

Vtip  (element  6  );  mass  8.32  mm/usec  ;  .087  g 

V  (element  100);  mass  .28  ;  2.78 

Jet  mass  66.3  g 

Jet  energy  .380  megajoules 

As  expected,  the  total  jet  energy  is  down.  The  tip  and  tai>  velocities 
are  nearly  that  of  the  baseline  charge,  as  well  as  the  tail  mass.  Total  jet 
mass  is  nearly  the  same,  but  tip  mass  is  about  50%  of  the  baseline  charge. 


Summary 


Charge 

Vtip 

Mtip 

Vtail 
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Mjet 

Ejet 

mm/usec 
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mm/usec 
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MJ 
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8.31 
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.33 

2.77 

71.6 
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5.67 
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.19 

2.76 
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RILC 

8.39 

.069 

.19 
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66.7 

.216 

SRILC 

8.34 

.070 

4.66 

.653 

32.1 

.616 

SRITLC 

8.45 
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36.3 
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Figure  6. 
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Conclusions 


In  spite  of  the  gross  assumptions  of  the  model,  the  results  appear  to 
be  intuitively  correct.  The  model,  perhaps,  does  not  give  accurate  absolute 
results,  but  the  relative  changes  in  the  jet  parameters  between  charges 
makes  sense.  In  this  limited  capacity,  the  model  has  some  utility  in  this 
analysis.  It  appears  from  the  results  that  of  the  six  jet  parameters 
evaluated  (jet  tip  velocity  and  mass,  tail  velocity  and  mass,  jet  mass,  and 
jet  energy)  at  most  four  parameters  can  be  reasonably  maintained  at  the 
expense  of  the  other  two.  However,  under  no  conditions  could  jet  tip  mass 
be  maintained  closely  to  that  of  the  baseline  charge.  It  is,  of  course,  up  to 
the  investigator  to  determine  which  parameters  are  critical  to  the 
experiment. 

The  best  results  appear  with  the  use  of  ring  initiation  and  a  tapered 
charge  This  may  cause  some  headache  in  manufacturing  and  assembly, 
but  can  certainly  be  done  with  liquid  explosives.  One  advantage  of  point 
initiation,  however,  is  that  more  head  height  takes  out  alignment  errors. 
Obviously,  more  precision  will  be  needed  for  ring  initiation. 
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Appendix  A 

Computer  Model  Output 


Legend: 

E#  Element  number 

Z  Element  axial  position  (mm) 

TJ  Element  time  to  reach  axis  (usee  after  charge  initiation) 

SM  Element  slug  mass  (grams) 

VS  Element  slug  speed  (mm/usec) 

JM  Element  jet  mass  (grams) 

VJ  Element  jet  speed  (mm/usec) 

BETA  Element  collapse  angle  into  stagnation  point  (degrees) 

VF  Element  flow  speed  into  stagnation  point  (mm/usec) 

FLAG  Integer  flag  if  flow  speed  exceeds  critical  velocity 
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